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Abstract 

This work addresses the response functions for infinite beams and plates with a force 
excitation at the origin of the coordinate system and its application to time reversal 
and time frequency analysis. The response function for Euler-Bernoulli beams is 
derived and for plates revisited. Interpretation concerning energy conservation and 
mobility are given. An possible alternative method to measure the coupling loss 
factors in SEA is sketched. 

The impulse response of a finite beam is measured in an experiment and compared 
with that predicted theoretically. The experimental data before the first reflection 
of the pulse showed good agreement with the theory. 

The response function is used to simulate a time reversed pulse. This numerical 
simulation is verified with an experiment. 
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1 Introduction 



This work addresses the response function for an infinite beam and its appli- 
cation to the time reversal in beams. The response function is an analytical 
function providing the responding velocity to ideally an impulse, but also to 
general excitation. 

The time reversal technique [1,2] is originally invented in the context of ultra- 
sonic sound and used mainly experimentally. The basic idea is first to record 
the impulse response of a source at many positions and then to playback this 
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response. The result is a strong focused signal at the source position. The only- 
prerequisite is reciprocity [3] . A theoretical study of the time reversal process 
in solids is in [4] and experimentally in [5]. The method is extended to other 
applications, e.g. damage detection [6]. 

Here the technique is applied to low frequency flexural waves in beams, for 
which the Euler - Bernoulli bending theory is valid. Whereas in [4] a general 
theoretical formulation is found, here the process is studied theoretically by 
means of the response function of a homogeneous beam. The advantage of this 
restriction is that an analytical solution can be studied. This gives additional 
insight in the phenomena and opens new possibilities for this specialised case, 
since it is not necessary to record the impulse. 

The application is restricted to beams, but nevertheless the corresponding 
function for an infinite plate is given. This shows that the concept can be 
easily extended to plates. 

The time reversal is mainly an experimentell technique. A theoretical approach 
that is similar to the presented can be found in [7,8] for the Timoshenko beam. 
In the study [7] the Green's function is obtained numerically whereas here an 
analytical function is studied. 

The response function itself is also useful for the analysis of the time - fre- 
quency distribution of a flexural impulse response. This application is the 
subject of a related publication [9] and is shortly introduced in section 4. 
The function itself was first derived for teaching purpose, since it is a instruc- 
tive example for dispersive wave propagation. In this context the functions are 
discussed with respect to energy considerations and the mobility in section 5. 
This provides interesting insight since well known phenomena are highlighted 
from a different point of view. 

The expectation is that given its basic character the response function should 
exist. An intensive search in mainly old scientific documents revealed the 
Green's functions for an initial deflection and velocity that is stated in the 
next section. From a phenomenological point of view there is maybe nothing 
to add since these functions are very similar. For a practical application there 
is a fundamental difference since for the later derived impulse response it is not 
necessary to know the initial deflection and velocity of the beam that result 
from the impulse. 



1.1 Green's function 

Green's functions for beams and plates were first derived by Boussinesq [10]. 
The governing equation for the deflection £ (x, t) of a beam is 
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where x is the distance from the source, ( = yB/m', where B is the flexural 
rigidity and m! mass per unit length. 

The Green's function is obtained by a Laplace transform of the differential 
equation (1) and its initial and boundary conditions. The Laplace transform 
and its application to transient vibration is well explained for example by 
Thomson [11]. With the initial conditions 



This function and its derivation can be found in several textbooks e.g. [12,13]. 
For applying the Laplace transform it is necessary to know the initial condi- 
tions. It is not obvious how to incorporate the excitation. For a plate, Boussi- 
nesq assumed that the initial deflection from a concentrating force acting at 
the origin of the coordinate system is £(r, 0) = This does not apply for 
the infinite beam as is shown later, equations see (14) and (15). 



2 Response Functions 

In this section the response functions for beams and plates are derived and 
discussed. In the following equation (1) is solved for the velocity v and not 
anymore for the deflection. The simple solution (13) that is derived in the 
following is only obtained if the velocity is chosen as the unknown variable. 
The approach is identical to that used in [14]. The difference to the example 
that is given in [14] is that the mobility is not altered by a sticking mass that 
impacts the beam. 

2.1 Semi-Infinite and Infinite Beams 

Consider a theoretical setup of a semi-infinite (x e [0, oo]) or infinite beam, 
which is excited at x = with an arbitrarily short force pulse F a . The pulse 
is modelled by the Dirac 5-function F a (t) = F S(t) whose Fourier transform 
is given by 



£(0,t)=V(t); £M) = 0; £(x,0) = 



(2) 



for the deflection £(x,t), it holds that 
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where H(T) is the Heaviside-Function. To find the response function to the 
impulse one may proceed with the boundary conditions for this elementary 
problem. The equations for the angular velocity w, bending moment M, shear 
force F y and velocity v of a beam which can be modelled with the Euler beam 
theory, equation (1) are 

A dv ~ B dw ~ dM . dF y 
w = — ) M — — — , F y = jum'v = -— ^. 5 

OX ]U) ox ox ox 

Further is the bending wave number fc 6 = u/cb, the bending wave velocity q, = 



^uj 2 B/m', the flexural rigidity B = EI y , E the elastic or Young's modulus, 
I y the geometrical moment of inertia. 

2.1.1 Boundary Conditions 

For a semi-infinite beam the waves propagate away from the excitation point 
which leads with v = Re{ve juJt } to 

v = v +e - jkiX + v +j e- kbX . (6) 

Herein the term v + is the amplitude of the far-field wave and v + j that of the 
near-field wave. At the free end the bending moment and shear force must 
vanish F(x = 0) = F a and M(x = 0) = 0. It follows with l+j = ^/2j that 



v(x, u) = " l ; f e~ jkbX + e~ kbX ) . (7) 
a/2j m'c b v ' 

In case of an infinite beam symmetry implies that 

v = v + e~ jklxl +v +j e~ klxl . (8) 
At x = the angular velocity is w(x = 0) = and F a /2 = F y (x = 0), therefore 

Hx,u) = -^^(e-^+e- k ^). (9) 



2m'cb\/2] 

This is the same result as in equation (7), except for the factor 1/4. 



2.1.2 Inverse Fourier Transform 

To obtain the response function v(x, t) one may use the real part of the inverse 
Fourier transform of equation (7) 

v(x,t) =Rej J v(x,u)e iuJt du\ , (10) 
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which with uo/yuj 2 = oyM ensures that the waves always propagate away 
from the excitation point. Upon substitution the expression reads 



v(x, t) - 

[ oo 

Re / 



F H(t) 
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where ( = ^ B/m' . This integral is solved with Mathematica's Integrate [] 
command with use of the hyper-geometric function p F q for t > 0. The hyper- 
geometric functions were found by Pochhammer[15]. Surprisingly he also solved 
the problem of bending waves in a circular beam without assuming that the 
radial movements can be neglected [16,17]. 

He gives also an exact expression for the bending wave velocity in this gen- 
eral case [18]. An overview of Pochhammer's work and its application towards 
transient response can be found in Miklowitz [19]. 

In the case of t < the integral vanishes due to the Heaviside-function. One 
thus finds that 



/ U — \f I LO I X U — \ 
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The pulse is an even function, so just the real part remains 




(13) 



It should be mentioned that to produce this result it is vital to use the near- 
and far-field terms in equation (11). For the origin at x — 0, it follows that 



v(x = 0,t) 



F H(t) 



(14) 



Since lim^^o ^ = 0, the cosine is unity and it follows for the velocity that 
\im x j^ov(x,t) — > oo. This relation can be found approximately in [14] for 
the example of a sticking mass. This example is maybe an indication that at 
least in Cremer's time the function (13) was not known, since a sticking mass 
is rather unusual and the force excitation has the advantage that it can be 
studied with an analytical function. 

The primitive function of equation (13) gives the deflection 
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where S(x) = Jq sin(7rt 2 /2) dt is the Fresnel-Sine-Integral. 



2.1.2.1 General excitation The Fourier transform if?(u) of a function 

4>(t) = 7 (t) x 4>{t) is 

q>(u) = —t(u) * (16) 

where * denotes the convolution. In case of an inverted Fourier transform there 
is no factor 1 / (2n) and it follows that 

„(*,*) = F.( t ).i ?1 gc«g. (17) 

For an arbitrary force F a (t), Im |_F a (c<j)| = does not hold. But F a (t) is always 
real and v(x,t) is real such that only the real part of equation (12) needs to 
be taken into account. It is a straightforward test to use the Dirac 5-function 
as the force and to obtain (13). 



2.2 Infinite Plate 



The derivation of the response function for the infinite plate is given by 
Crighton [20,21]. Nevertheless this problem is revisited, since the equation 
is used in the following. The starting point is the propagation function of an 
infinite plate 

v(u>,r) = ^=L (H { i\k b r) - H^\-jk b r)) , (18) 

with the Hankel function of the second kind Hq 2 \x) = Jq(x) — jY (x), the 
flexural rigidity of a plate B' = 12 fi_ u \ > v the Poisson ratio and m" the mass per 
unit area. The inverse Fourier transform can be simplified with the following 
identies of the Hankel function. For x > are 



H o 2 \-jx) = -jRe{Y (-jx)} , J (x) = J (-x), 
Im{y (-^)} = 2J (-x), Im{J } = 0, 

Re{y (^)} = Re{Y (-jx)}, Im{y ^)} = 0, 
Re{Y (x)} = Re{Y (-x)} 



(19) 



With (19) the real part of the inverse Fourier transform of equation (18) 
reduces to 

v(r, t) = ~^=jj J U^Rr) cosM) du. (20) 
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Hence, the response function for the infinite plate is 

F . r 2 



v(r, t) = ^ sin — . (21) 



3 Dispersion factor and number 



An important difference between the bending wave and, for example, a lon- 
gitudinal wave is that the group velocity is frequency dependent and so the 
response to an impulse is spreading into the different spectral fractions of the 
pulse. The term 

* = ^ (22) 

is the factor that controls this spreading and is called dispersion factor. 
Whereas the dispersion factor is a time value the nondimensional term 

Di = X U^^ (23) 

is called dispersion number. In an experiment that should reveal the influence 
of the dispersion it is necessary to choose a structural sample with a high 
dispersion number. The maximum frequency is impotant since the higher the 
frequency the more the impulse will spread. 

A high dispersion is the reason for choosing a thin plate and a slender beam. 



4 Time - Frequency Analysis 



The time-frequency analysis of dispersive waves is widely studied and can be 
found in several publications. This section is only meant to give an overview 
of applications of the response function in this area. 
One may insert x = c g t in equation (13), to obtain 

, x FqH(t) PT . 

v(u,t) = y - * —cosut, (24) 

where the bending wave group velocity is c g = 2q, = 2\^u(. This shows 
the known fact that the frequency content of the bending wave is travelling 
with its particular group velocity. It is the theoretical basis for the time- 
frequency analysis of dispersive waves. For example Kishimoto et. al. [22] show 
the validity of their approach with two neighbouring frequency components. 
Often the wavelet transform [22,23] is used to extract the frequency dependent 
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arrival time t a = c g /x. 

The time and frequency dependent energy is obtained from the signal by 
means of the wavelet transform. The proceduce to extract the arrival time is 
to find the time value with the maximum energy at a fixed frequency. It is 
not possible to find the frequency value with the maximum energy at a fixed 
time, since this is e.g. affected by the non-constant amplitude F of a real 
application. 

An extention of the described method of time- frequency analysis is done by 
Vries et. al.[24]. The basic idea is to remove the dispersion by applying the 
response function. The response function itself is obtained by a numerical 
inverse Fourier transform. 

Another approach is to use the analytical response function (13) as a wavelet 
mother function. The dispersion factor is then obtained by the scaling factor 
of the wavelet with the highest energy. This approach can be found in [9] and 
is applied to the signals obtained in the experiments that are presented in the 
following. 

That the approach above is valied can also be shown with the instaneous 
frequency cu(t) of an almost periodic functions, see Bochner [25]. Bochner [26] 
also showed that these functions are solution to the wave equation and also 
for equation (1) that is used here. It holds that 

cos<p(t)^u(t) = ip'(t) = -^. (25) 



5 Energy Conservation and Mobility 



The term y 7 in equation (13) stems from the general energy conservation 
scaling of a function f(x), which is 




One may interpret equation (24) so, that the pulse while traveling along the 
beam is scaled by the travel time such that the energy of the pulse is conserved. 
In the case of a plate the corresponding term is |, which is due to the fact 
that the pulse propagates cylindrically and not as a plane wave. The radius of 
the cylinder follows ar~\/t dependence. The complex mobility is defined by 

Y = v/F a . (27) 

If the velocity and the force are not at the same position it is called transfer 
mobility. 

For a semi-infinite beam the transfer mobility to a position in the far-field 
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(ki,x 3> 1) is obtained as 



Y(u) = = , -e- ]kbX . (28) 

1 ; v^Cm'(l+j) 1 ; 



The decrease in magnitude over frequency is given by 

w - hM>- (29) 

That this corresponds to the solution in the time domain is discussed in the 
following. Consider the envelope of equation (13) to be 



FpHjt) 2 

Venv = ( 30 ) 

m y nQt 

The Fourier transform of equation (30) produces with FT{^J\Jt} = \Jir/u> 
exactly the same mobility as equation (29). This indicates that the remaining 
cosine-term in equation (13) is not affecting the magnitude of the mobility in 
the far- field. From equation (33) follows that the far- field condition (k^x ^> 1) 
for equation (13) is di/t ^> 1. With this prerequisite it follows that the envelope 
(30) defines the amplitude of a frequency component given by the remaining 
cosine-term, like in equation (24). For low frequencies di/t <C 1 it is not 
possible to analyse the two terms independently, since the prerequisite for 
almost periodic functions is violated [25]. Nevertheless the Fourier transform 
of equation (13) results obviously in equation (7). 
With this interpretation one can define a mobility in the time domain 

Yt = v(t)/F a (u = . (31) 

This definition is followed in the experimental section 6. It also leads to a 
possible application that is discussed in the following. 



5.1 Application to SEA 



It is shown that the mobility is trivially related to the response function. This 
opens a possible application in the context of the Statistical Energy Analysis 
SEA [27]. To measure the coupling loss factors used in a SEA model a stan- 
dard method is to excite each subsystem in turn and measure the response of 
all the subsystems for each source [27,28]. 

The elaborated interpretation opens an alternative method: excite a subsystem 
A with an impulse and measure the transient impulse reponse of the neigh- 
bouring subsystem B. The transmission coefficient can be easily obtained with 
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the energy of the impulse and its response. The only prerequisite is that the di- 
rect transmission can be clearly distinguished from the early reflections. That 
this is in principle possible is shown in the experiments in section 6, where the 
response function for an infinite beam is verified with a finite beam. 



6 Experimental results 

Measurements are carried out on a thin acrylic plate and a slender aluminum 
beam for different distances and configurations. The results show the same 
tendency. For the sake of brevity just one typical measurement of the beam 
and the plate is presented and discussed. 



6. 1 Beam 



The dimensions of the beam are a diameter of d = 6mm, a length of I = 3m 
and a density of p = 2830kg /m 3 . The material parameters for aluminium 
are tested with resonance frequencies of the beam and a elastic modulus 
E = 68, OGPa is obtained. 

The beam is clamped at both ends. The velocity is measured with a laser- 
vibrometer with a sampling frequency of 50kHz placed at 0.5m from the end. 
The beam is excited by means of an impacting hammer equipped with a force 
transducer at a distance of x = 1.50m from the point of measurement. 
Since the real beam is not infinite only the first passage of the pulse is consid- 
ered. The beam velocity is plotted in figure 1. The time axis is started at the 
maximum of the force signal minus the delay of the laser- vibrometer of 1.2ms. 
One may recognise that already after 3ms the reflections from the clamped 
end are visible. 

The theoretical curve is calculated by means of equation (13) with a value of 
di = 0.079. This corresponds quite well with the value obtained with equa- 
tion (22) of d{ = 0.076. The actual value is extracted with a method that is 
discussed in a subsequent publication. The frequency range of the theoretical 
curve is f min = 124Hz to f max = hAkHz. The measured curve is corrected by 
means of the frequency distribution of the theoretical curve that is obtained 
from equation (33) and the power spectrum of the measured force impulse, 
like defined in equation (31). 

The curves are normalised with their maximum value, since in this context 
the distribution of amplitude and frequency over time is of interest, but not 
the absolute value. 
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Fig. 1. Impulse response of a beam, normalised velocity 
v(x = 1.5m, t)/\v(x = 1.5m, t)\ max : measured (plus), theoretical (diamond) 
with equation (13), measured and corrected with the power spectrum of the 
impulse (circle) 

6.2 Plate 

The dimensions of the plate are a thickness of d — 2mm, a length of I = 2.05m 
and a width of b = 1.52m. The material parameters provided by the manufac- 
turer are elastic modulus E = 3.3GPa, p = 1190kg/m 3 and a Poisson's ratio 
v = 0.37. 

The whole plate is suspended in a frame. The velocity is measured with a 
laser-vibrometer at positions at least 0.5m from the frame. The excitation 
point is at a distance of 0.5m from the response position. 
The velocity in figure 2 is obtained in almost the same manner as in the beam 
experiment. The theoretical curve is calculated with equation (21) and a value 
of di = 0.0544 in a frequency range from f min = 180Hz to f max = 5.2kHz. 



7 Time reversal 

To generate a defined pulse at a certain distance from the point excitation one 
may develop the time reversal from equation (13). As an example, a Dirac like 
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Fig. 2. Impulse response of a plate, normalised velocity 
v(r = 0.5m,t)/\v(r = 0.5m,t)\ max '- measured (plus), theoretical (diamond) 
with equation (21), measured and corrected with the power spectrum of the 
impulse (circle) 

pulse is requested at a position a; on a structure. By considering the force 



Fa(t) 



"s/t-max t 





cos 



4C(t« 



-tr- 



ior t r 



< t < t r , 



otherwise . 



(32) 



a pulse resembling a ^-function can be realized. The generated pulse will not 
be a perfect 5-function that consists of the whole frequency spectrum, but will 
be a band-filtered version. 

With equation (25) the frequency range of equation (32) is given by 



fmax/r, 



X 



(33) 



mini max 



In a numerical experiment the velocity is calculated with equation (17) and 
shown in figure 3. 



7. 1 Experimental Results 



The studied beam is the same like in section 6.1, but now the beam is clamped 
at one end and hanging on a twine on the other. The measured velocity is anal- 
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t 

Fig. 3. Velocity of a beam excited with a time reversed pulse defined in equation 
(32), £ = 1, x = 3. The darkness is given by \v\. 

ysed in the same manner as in the previous experiments. 
The generation of the desired force excitation is done with a magnetic force 
transducer. The dispersion factor for the reversed pulse (32) is chosen to 
di = 0.0054 with the given material properties of the beam the location of 
the pulse is at Xi mp = 0.4m. The velocity is measured with a laser vibrometer 
in a distance of x meas = 1.13m. It is assumed that from the receiver position 
a virtual source at a distance of x v = 0.73m is seen. 

The measured velocity is analysed with the method that is discussed subse- 
quent publication and a dispersion factor of di = 0.0181 is extracted within 
a frequency range between f m ax = 7.8kHz and j m in = 1.0kHz. The value of 
di = 0.0181 corresponds with the distance of the virtual source of x v = 0.73m. 
The normalised velocity is plotted in figure 4. The theorectical curve is cal- 
culated with the dispersion factor of di = 0.0181. The results are not that 
clear as in the previous experiments. It is assumed that this is due to the 
comparable weak magnetic force transducer. 



8 Concluding Remarks 

The main focus of this study is on the derivation of the impulse response (13) 
and its application. It is a simple and instructive example of dispersive wave 
propagation. 
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Fig. 4. Normalised velocity of a beam excited with a time reversed pulse defined in 
equation (32) v{x = 1.13m, t)/\v{x = 1.13m, t)\ max : measured (dashed), theoretical 
(solid) with equation (13) for a distance of x v = 0.73m 

The applications of the impulse response function for time reversal is shown 
in theory and experiment. The attempt differs since an analytical function is 
used and not the measurement impulse response. 

Other applications in the fields of time frequency analysis and SEA are in- 
dicated. This applications are rather roughly sketched and are subjects of 
ongoing research. 

Another natural extension is to study the conversion from transient to sta- 
tionary motion. 
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